Influence of electron correlations on ground-state properties of III-V semiconductors 
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Lattice constants and bulk moduli of eleven cubic III-V semiconductors are calculated using an 
ab initio scheme. Correlation contributions of the valence electrons, in particular, are determined 
using increments for localized bonds and for pairs and triples of such bonds; individual increments, 
in turn, are evaluated using the coupled cluster approach with single and double excitations. Core- 
valence correlation is taken into account by means of a core polarization potential. Combining the 
results at the correlated level with corresponding Hartree-Fock data, we obtain lattice constants 
which agree with experiment within an average error of —0.2%; bulk moduli are accurate to +4%. 
We discuss in detail the influence of the various correlation contributions on lattice constants and 
bulk moduli. 
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I. INTRODUCTION 



II. HARTREE-FOCK CALCULATION 
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Ground state properties of solids are most often calcu- 
lated within the framework of density-functional theory 
(DFT), where the exchange and correlation contributions 
are determined within the local-density approximation 
(LDA)El. These methods yield good results for lattice 
constants and bulk moduli but do not provide for many- 
body wave functions. Both, the influence of exchange and 
correlation are covered in an implicit way, thus, a system- 
atic improvement towards the exact results appears to be 
difficult. 

Hartree-Epck self-consistent-field (HF-SCF) calculations 
for solidscl, while lacking electron correlations as per its 
definition, have the merit of treating non-local exchange 
exactly. They supply a good starting point for subse- 
quent correlation calculations provided correlations are 
not too strong. Infinite systems require the use of size- 
consistent approximations for electron correlations. In 
Ref. ^ an approach based on cumulants was discussed 
which ensures size-consistency at any level of approxi- 
mation. We use it here in order to work with local oper- 
ators, i.e., local one- and two-particle excitations. These 
are best suited for describing the local correlation hole 
which the electrons carry with them. The matrix el- 
ements which appear in the calculation of the ground- 
state energy are evaluated using clusters instead of solids 
thereby applying an incremental methodu. The same 
method has previously been used to calculate_ground- 
statc properties of elementary semiconductorsU|-as well 
as the cohesive energy of cubic III-V compoundsO. 
In the present paper we apply this method to eleven III-V 
compounds in zinc-blende structure and investigate the 
influence of electronic correlations on lattice constants 
and bulk moduli. In addition to the valence correlations 
the core polarization, which roay have large influence on 
the bond lengths of moleculesa, is taken into account for 
the solid in a systematic way. 



Before discussing the influence of electron correlation 
effects in solids, reliable HF-SCF calculations are a neces- 
sary prerequisite. We performed HF ground-state calcu- 
lations for the IH-V compounds using the program pack- 
age Crystal92lI. The same basis sets and pseudopoten- 
tials as in Ref. || are usedB. 

The HF lattice constants have been determined from the 
minima of the HF ground-state energy curves in the fol- 
lowing way. We varied the lattice constants in steps of 
1% of the experimental values and calculated six points of 
the potential curves in each case. The minima were eval- 
uated from a quadratic fit to these points. A narrower 
spacing and higher order fits were tested but yielded only 
a marginal change in the position of the minimum and 
the curvature. The results for the HF lattice constants 
are listed in Table I. 

The bulk moduli, which describe the response of the solid 
to a homogeneous pressure, have been evaluated at the 
experimental lattice Constantsa. The HF values for the 
bulk moduli are listed in Table II. 



III. CORRELATION EFFECTS 



A. Core polarization 



The description of core electrons by means of pseu- 
dopotentials implies a neglect of static and dynamic core 
polarization. The first part is an SCF effect, while the 
latter part is related to core-valence correlation. Espe- 
cially if closed d shells are described by pseudopotentials, 
the influence of the core polarization on lattice constants 
and bulk moduli is significant although the influence on 
cohesive energies is negligible. The computational effort 
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for an explicit treatment of the core electrons and their 
correlations would be very high; thus, we simulate these 
effects in our calculations by means of a core polarization 
potential (CPPp. ft describes the dipole interaction be- 
tween the core and the valence electrons: 



Vcpp = -£t^a/a; 



(1) 



here A numbers different cores, a\ is the corresponding 
dipole polarizability, and f\ is the field at the site of 
the core generated by valence electrons and surrounding 
cores: 
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f\i is the distance between valence electron i and core 
A, Ny is the number of valence electrons. The cut-off 
parameter 5\ is necessary in order to remove the singu- 
larity of the dipole interaction at r\i = 0. We took the 
parameters a\ and d\ from Ref. [l(], where the CPP was 
adjusted in atomic calculations to the spectra of single- 
valence-electron ions. 

In order to obtain information on the influence of CPP's 
on solid-state properties, we choose a finite fragment of 
the zinc-blende structure - a A4I4H18 cluster - where the 
dangling bonds are saturated with hydrogen atoms. The 
X — H and the Y — H distances, respectively, n are opti- 
mized in CCSD calculations for AVHg clusters!!. For hy- 
drogen we choose Dunning's double-zeta basisllil without 
the p polarization function. For the other elements, large- 
core pseudopotentials in conjunction with corresponding 
(4s4plG?)/[3s3pfd] valence basis sets are usecOfi. 
We first want to check the local character of the CPP. 
Test calculations for GaAs are performed in the follow- 
ing way. The CPP for Ga is only provided at the in- 
nermost Ga atom and the energy gain due to the CPP, 
A£cpp(Ga), is calculated at the_.HF-SCF level with 
the program package Molpro94E3. In the next step 
we provide the CPP only at the innermost As atom, 
and Ai?cpp(As) is determined. Providing the CPP's 
at the inner Ga and the inner As atoms simultaneously, 
we finally calculate A£"cpp (Ga, As) , and we find that 
the non-additive part Ai5cpp(Ga, As) — (AEcpp(As) + 
AEcpp(Ga)) is negligibly small. 

With this knowledge in mind, we determine an approx- 
imation for the core-polarization energy of the solid per 
unit cell by calculating AEcpp(X,Y) in the A4I4H18 
cluster as described above. Varying all X — Y bonds 
in the X4Y4H1S cluster, we obtain the dependence of 
AEcpp(X,Y) on the lattice constant a. The results at 
this level are listed in Table I for the lattice constants and 
in Table 11 for the bulk moduli; they include core- valence 
correlation (through the use of the CPP) but do not ac- 
count for the coupling between core-valence and valence 
correlation, which is dealt with in the next subsection. 



B. Valence correlations 



As pointed out in the introduction, the matrix ele- 
ments which appear in the calculations of the ground- 
state energy are obtained from replacing the solid by suf- 
ficiently large clusters. With these matrix elements the 
correlation contribution to the ground-stata-energy can 
be partitioned into a sequence of incrementsaa which are 
successively evaluated. The procedure and the computa- 
tional details for the III-V compounds are described in 
Ref. In order to account for the valence/core- valence 
correlation-energy coupling, we included CPP's on every 
X- and V-atom and used the HF-SCF ground state with 
CPP as reference state for the correlation calculation. 
For the evaluation of the largest increments (one-bond 
and up to next-nearest-neighbour two-bond) the basis set 
was extended to (4s4p2<il/)[I Applying this method we 
determined the correlation energies as function of the lat- 
tice constant. The calculated lattice constants and bulk 
moduli are listed in Table I and Table II. 



IV. RESULTS AND DISCUSSION 



We have calculated the lattice constants and the bulk 
moduli for eleven cubic III-V semiconductors, at differ- 
ent theoretical levels. A HF treatment overestimates the 
lattice constants by up to 2% , with the errors increasing 
for systems with heavier atoms. Apparently, it is due to 
the lack of electron correlations which describe the in- 
stantaneous motion of the electrons, that the structural 
minimum is found at larger volumes as compared to ex- 
periment. This is opposite to what is usually found for 
molecules composed of first- and second-row atoms in 
quantum-chemical calculations. There, the lack of non- 
dynamical correlation (left-right correlation), yields an 
artificial upwards shift of the potential-energy curve for 
large R and consequently leads to too short SCF bond 
lengths. Of course, this effect is also present in the solids 
considered in this paper, but apparently is overridden by 
dynamical correlation effects (cf. below). However, quite 
in common with quantum-chemical experience for molec- 
ular force constants, the bulk moduli are overestimated 
by up to 30%, at the HF level. The materials seem to be 
harder (large bulk modulus) because the instantaneous 
answer of the electrons to pressure is neglected. 
In the next step we take into account the core polariza- 
tion. As described in Section III. A, the core polarization 
is neglected when using pseudopotentials. Pictorially 
speaking, the instantaneous deformation of the spheri- 
cal cores in the field of the valence electrons is omitted. 
Modelling this effect with a core polarization potential, 
we achieve a significant reduction of the lattice constants 
(by up to -2%). The calculated values for the bulk mod- 
uli decrease by up to 15% depending on the materials 
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(heavier atoms are easier to polarize). 
The overall effect of the valence correlations determined 
with the method of local increments leads to a reduction 
of about 1% in the lattice constant and about 10% in the 
bulk modulus. 

To get a deeper understanding of the nature of the in- 
terplay of the various correlation effects, we performed 
calculations for GaAs with different basis sets (see Ta- 
ble III). The interatomic (non-dynamical) part of the va- 
lence correlation energy is determined using a minimal 
(4sAp)/[lslp] basis set. It enlarges the lattice constant by 
more than 2%, because in the many-body wave function 
only excitations into anti-bonding orbitals are possible. 
Intra-atomic (dynamical) correlations (calculated with a 
larger than minimal basis) decrease the lattice constant 
because more low lying excitations are mixed into the in- 
teratomic many-body wave function within the enlarged 
configuration space. The electrons can avoid each other 
much better. Increasing the basis further (to 4s4p2dl/)El, 
the intraatomic correlations are increased and lattice con- 
stants and bulk moduli are decreased. 
A comparison between valence-correlation effects with 
and without CPP has also been made for GaAs (see Ta- 
ble III). This way, the coupling between valence and core- 
valence correlation can be investigated. While the inter- 
atomic part of the valence-correlation energy is virtually 
unchanged when including the CPP, the intraatomic one 
is reduced in the calculation with CPP. The reason for 
this effect, however, is not so much a reduction of dynam- 
ical valence correlation but rather a reduction of dynam- 
ical core-polarization when angular correlation around 
the atoms is introduced into the many-body wave func- 
tion: angular correlation leads to a preference for valence 
electrons occupying opposite positions with respect to an 
atomic core so that their polarizing effect is reduced as 
compared to an independent-particle (uncorrelated) de- 
scription. The major part of this effect comes from s, 
p and d functions, the additional bond- length reduction 
due to / functions is nearly the same in the calculation 
with and without CPP. 

For the lattice constants we reach a very good agreement 
with experiment. All calculated results are slightly too 
small (by -0.2% on the average), but they are within the 
error bars given by the pseudopotentials and the CPP's 
(~ 0.03 A). We have checked the quality of the chosen 
basis set for the Ga — As bond length in GaAsHg . There, 
we see nearly no change (<0.1%) when extending the ba- 
sis to an even larger one (6s6p3d2flg). The average error 
of the calculated bulk moduli is about +4% . 
For comparison, we have also listed in Table I and Table 
II results from the literature which have been obtained 
with LDA. Lattice constants are underestimated within 
the LDA treatment by up to 2.8%. Bulk moduli mostly 
agree well with experiment. 



V. CONCLUSION 

We have determined the lattice constants and the bulk 
moduli of eleven cubic III-V semiconductors. At the HF- 
SCF level, the lattice constants are overestimated by up 
to 2%, the bulk moduli by up to 30%. A significant 
reduction towards the experimental values is obtained 
when applying core polarization potentials. Valence cor- 
relation has been calculated at the CCSD level, using the 
method of local increments. For the lattice constants two 
opposite trends are seen: interatomic correlations enlarge 
the lattice constant, while intraatomic ones decrease it. 
The overall trend is a reduction of the lattice constants, 
and the final results agree very well with experiment (av- 
erage error ta-0.2%). Valence correlations substantially 
decrease the bulk moduli, and again the agreement with 
experiment is remarkable good (average error «+4%). 
Due to the availability of many-body wave functions in 
our scheme, the influence of correlations on the lattice 
constants and the bulk moduli can be discussed in de- 
tail. 
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TABLE I. Lattice constants (in A ) at different theoretical levels (jgf. text); deviations from experimental values are given in 
parentheses; the experimental values measured at room temperaturefj have been extrapolated to zero Kelvin. 



ClHF 



flHF+CPP 



tSHP+CPP+corr 



IlLDA 



4.5836 
4.8276 
5.5348 
5.7405 
6.2606 
5.5298 
5.7546 
6.2120 
5.9527 
6.1524 
6.5925 



(+ 1.0%) 
(+1.1%) 
(+1.3%) 
(+1.5%) 
(+2.1%) 
(+1.5%) 
(+1.9%) 
(+2.2%) 
(+1.5%) 
(+1.6%) 
(+1.9%) 



4.5679 
4.8011 
5.5016 
5.6944 
6.1954 
5.4636 
5.6730 
6.1034 
5.8559 
6.0521 
6.4717 



(+0.7%) 
(+0.5%) 
(+0.7%) 
(+0.6%) 
(+1.0%) 
(+0.3%) 
(+0.4%) 
(+0.4%) 
(-0.2%) 
(±0.0%) 
(±0.0%) 



4.5322 
4.7644 
5.4398 
5.6428 
6.1307 
5.4274 
5.6472 
6.0740 
5.8296 
6.0437 
6.4615 



(-0.1%) 
(-0.3%) 
(-0.5%) 
(-0.3%) 
(±0.0%) 
(-0.3%) 
(±0.0%) 
(-0.1%) 
(-0.6%) 
(-0.2%) 
(-0.2%) 



4.777 
5.421 
5.62ai2 



5.358 

5 ' 939 T3 
5.7021 LS 

5.9019 181 



(-1.4%) 
(± 0.0%) 
(-0.8%) 
(-0.7%) 

(-1.7%) 
(-2.5%) 
(-2.3%) 
(-2.8%) 
(-2.5%) 
(-2.0%) 



TABLE II. Bulk moduli (in Mbar) at different theoretical levels (cf. text); deviations from experimental valueaj are given in 
parentheses. 





Bhf 


Bhf+ 


CPP 


Bhf+ 


CPP+corr 


BP 


1.95 


(±12%) 


1.92 


(±11%) 


1.87 


(+8%) 


BAs 


1.62 




1.53 




1.51 




A1P 


1.01 


(±17%) 


0.98 


(±14%) 


0.92 


(+7%) 


AlAs 


0.90 


(±17%) 


0.86 


(±12%) 


0.82 


(+7%) 


AlSb 


0.65 


(±12%) 


0.61 


(±5%) 


0.57 


(-1%) 


GaP 


1.10 


(±21%) 


0.98 


(±8%) 


0.93 


(+2%) 


GaAs 


0.87 


(±13%) 


0.82 


(±6%) 


0.77 


(0%) 


GaSb 


0.63 


(±12%) 


0.63 


(±12%) 


0.59 


(+6%) 


InP 


0.79 


(±10%) 


0.71 


(-1%) 


0.66 


(-8%) 


InAs 


0.74 


(±28%) 


0.69 


(±19%) 


0.64 


(+10%) 


InSb 


0.58 


(±26%) 


0.52 


(±11%) 


0.48 


(+4%) 



n 



1.7312 



0.978a 
0.77F 
0.799 
0.761^ 
0.619C 
0.4774L! 



(± 0%) 



(±7%) 
(± 0%) 
(±43%) 
(±6%) 
(±7%) 
(±4%) 



TABLE III. Influence of correlations on lattice constants (in A) of GaAs, employing different basis sets (cf. text); deviations 
from experimental values are given in parentheses. 





HF 


CCSD, min. basis 


CCSD,3s3pld 


CCSD, 4s4p2dl/ 


a (with CPP) 


5.6730 (±0.4%) 


5.7910 (±2.5%) 


5.7084 (±1.1%) 


5.6472 (±0.0%) 


a(without CPP) 


5.7546 (±1.9%) 


5.8810 (±4.1%) 


5.7645 (±2.1%) 


5.7015 (±0.9%) 
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